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ABSTRACT 


This  report  deals  with  the  prediction  of  service  volumes  (i.  e., 
service  ranges  in  three  dimensions)  for  VOR  and  TACAN  air  naviga¬ 
tion  facilities  in  the  presence  of  noise  and  interfering  stations.  Propa¬ 
gation  mechanisms  applicable  to  the  frequency  ranges  employed  are 
discussed  first,  together  with  calculation  of  transmission  loss  and  its 
variability.  Second,  the  statistical  treatment  of  the  interference 
problem  is  explained.  Finally,  it  is  shown  how  service  volumes  are 
determined  based  on  transmission  loss  and  its  variations,  specified 
or  assumed  criteria  on  signal-to-noise  and  signal-to-interference 
ratio,  and  the  spatial  relations  of  the  aircraft,  the  desired  facility, 
and  the  undesired  facility.  The  results  of  the  entire  study  are  pre¬ 
sented  in  graphical  form. 

Detailed  procedures,  mathematical  formulas,  and  computer 
programs  used  are  discussed  in  the  Appendix. 


1.  INTRODUCTION 

Increasing  air  traffic  density  together  with  fast,  high-flying  jets 
have  made  the  use  of  reliable  air  navigation  aids  more  imperative  than 
ever  before.  In  expanding  the  present  complex  of  navigation  aids  to 
meet  future  demands  consideration  must  be  given  to  potential  inter¬ 
ference  between  facilities  operating  on  the  same  or  on  adjacent 
channels.  The  amount  of  interference  is  a  function  of  the  desired-to- 
unde sired  signal  ratio  at  the  aircraft  antenna  terminals;  as  both  signals 
vary  with  time  and  aircraft  location,  the  ratio  varies  as  well,  and  the 
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interference  becomes  dependent  on  time  and  location.  Due  to  the 
nature  of  radio  wave  propagation  in  the  frequency  ranges  used,  the 
variations  of  the  received  signals  and  of  the  interference  ratios  are 
best  described  statistically.  The  large  number  of  possible  conditions 
dictates  the  use  of  a  digital  computer  with  programs  that  take  into 
account  all  variables  as  well  as  the  fixed  equipment  parameters. 

The  navigation  aids  treated  here,  namely  VOR  and  TACAN 
operate  in  the  100  and  1,  000  Me/ s  frequency  ranges,  respectively. 

In  these  ranges,  propagation  of  radio  frequency  energy  is  affected 
principally  by  the  lower  atmosphere  (the  troposphere),  specifically 
by  variations  in  the  refractive  index  of  the  atmosphere.  The  terrain 
along  and  in  the  vicinity  of  the  great  circle  path  between  transmitter 
and  receiver  also  plays  an  important  part. 

Within  the  last  decade  a  number  of  methods  and  procedures  have 
been  developed  to  calculate  field  strength  and  its  variability.  In  this 
report  procedures  are  followed  which  have  been  used  at  the  National 
Bureau  of  Standards  to  predict  in  a  statistical  manner  the  effects  of 
terrain  and  atmosphere  on  the  variability  of  field  strength  and  on  the 
performance  of  radio  systems;  see  for  example  Barsis,  Norton,  Rice, 
and  Elder  [1961]  ,  and  NBS  Report  67  67  [1961]  .  It  is  also  convenient 
to  use  the  concept  of  transmission  loss  [Norton,  1953  and  1959],  which 
is  the  ratio  of  power  transmitted  to  the  power  available  at  the  receiv¬ 
ing  antenna  terminals,  usually  expressed  in  decibels.  Methods  used 
for  its  calculation  as  a  function  of  path  length,  terminal  heights,  and 
carrier  frequency  are  discussed  below. 

All  fixed  and  variable  parameters  for  either  the  VOR  or  the 
TACAN  system  are  assembled  into  essentially  one  general  computer 
program  which  yields  signal-to -noise  and  signal-to-interference  ratios 
for  given  probability-of-service  values  as  functions  of  aircraft  location 
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in  relation  to  the  desired  and  the  undesired  transmitting  facility  and 
equipment  parameter  s .  From  these  results,  boundaries  of  service 
volumes  have  been  established  corresponding  to  points  in  space  where 
the  specified  desired-to-undesired  signal  strength  ratios  will  be 
exceeded  with  a  probability  of  0.9  5. 

For  the  time  being,  only  one  interfering  station  has  been  con¬ 
sidered  in  any  particular  case,  and  the  study  is  restricted  to  conditions 
along  a  common  great  circle  path  between  the  desired  station,  the 
aircraft,  and  the  undesired  station. 

Z.  TRANSMISSION  LOSS  CALCULATIONS 

Figure  1  shows  a  typical  contiguration  of  the  aircraft  (represent¬ 
ing  the  receiving  terminal),  a  desired  navigational  transmitting  facility, 
and  an  undesired  navigational  transmitting  facility.  All  three  are 
aligned  along  a  great  circle  path,  and  for  simplicity  assumed  to  be 
above  a  smooth  surface.  In  the  example  drawn,  the  aircraft  is  within 
the  radio  horizon  of  the  desired  facility,  but  beyond  the  radio  horizon 
of  the  interfering  station.  The  distances  along  the  great  circle  path 
from  a  point  vertically  below  the  aircraft  to  the  desired  and  the  un¬ 
desired  station  are  denoted  by  d^  and  d  ,  respectively.  The  air¬ 
craft  is  at  a  height  h^  above  the  terrain.  The  angle  9  between  the 
horizon  rays  from  the  aircraft  and  the  interfering  station  is  an  impor¬ 
tant  parameter  in  the  calculation  of  transmission  loss  for  beyond-line- 
of-sight  paths  [Norton,  Rice,  and  Vogler,  1955]  .  For  the  assumption 
of  a  smooth  spherical  earth,  9  is  also  the  distance  between  radio 
horizons  expressed  in  angular  measure,  and  it  is  therefore  called  the 
"angular  distance,  "  as  indicated  on  Fig.  1.  As  is  customary  in  the 
analysis  of  tropospheric  propagation  paths,  first  order  allowance  for 
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the  refraction  effects  of  the  atmosphere  is  accounted  for  by  assuming 
an  earth  radius  4/3  times  the  actual  radius.  This  permits  radio  rays 
to  be  drawn  as  straight  lines,  and  simplifies  the  determination  of 
geometric  parameters. 

Transmission  loss  calculations  are  accomplished  in  three  steps 
for  the  application  to  air  navigation  problems,  as  follows: 

(a)  A  reference  value  of  basic  transmission  loss  is  calculated 
in  accordance  with  methods  given  in  NBS  Report  67  67 
[1961]  .  This  reference  value,  as  explained  by  Norton, 
et  al.  [1955],  represents  the  conditions  on  an  average 
winter  afternoon  in  the  northern  temperate  zone.  More 
specifically,  it  is  the  median  of  all  hourly  medians  within 
"Time  Block  2,  "  representing  the  hours  1:00  p.m.  to 
6:00  p.m.,  November  through  April.  The  term  "basic 
transmission"  loss  means  that  for  the  purpose  of  these 
calculations  the  antennas  at  the  terminals  are  considered 
to  be  isotropic  with  unity  power  or  voltage  gain  in  all 
directions.  Actual  antenna  characteristics  are  brought 
in  later  on  when  converting  from  basic  transmission  loss 
to  transmission  loss. 

Within  the  radio  horizon,  basic  transmission  loss  is 
calculated  using  geometric  optics  methods  including  inter¬ 
ference  between  the  direct  and  the  ground- reflected  ray. 
For  the  VOR  propagation  model  in  the  110  Mc/s  frequency 
range,  a  specular  reflection  was  assumed  with  effective 
reflection  coefficients  on  the  order  of  0.9.  In  the  case  of 
TACAN  with  frequencies  in  the  1100  Mc/s  range  it  is  more 
appropriate  to  assume  that  the  ground- reflected  ray  is 
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made  up  of  a  large  number  of  components  having  random 
relative  phase.  The  total  contribution  from  these  compon¬ 
ents  can  be  represented  statistically  by  a  Rayleigh  distri¬ 
bution  [Decker,  19  57;  McGavin  and  Maloney,  19  59]  .  Norton, 
Vogler,  Mansfield,  and  Short  [1955]  have  shown  how  this 
Rayleigh-distributed  reflected  component  can  be  combined 
with  a  fixed  component  representing  the  direct  ray  in  order 
to  arrive  at  the  resulting  distribution  of  basic  transmission 
loss  at  the  receiver. 

Beyond  the  radio  horizon,  basic  transmission  loss  is 
calculated  using  smooth-earth  diffraction  or  forward 
scatter  models,  depending  on  the  path  distance  involved. 

The  diffracted  field  decreases  very  rapidly  for  distances 
beyond  the  radio  horizon,  especially  at  the  TACAN  fre¬ 
quencies,  so  that  the  forward  scatter  model  is  the  more 
important  one.  Calculations  for  both  models,  and  the 
method  of  properly  combining  diffraction  and  scatter  fields 
if  they  are  of  comparable  magnitude,  are  based  on  pro¬ 
cedures  contained  in  NBS  Report  67  67  [1961]  . 

(b)  Long-term  variations  in  basic  transmission  loss  are  esti¬ 
mated  by  means  of  time  variability  functions  V(p,  9)  or 
V(p,  d)  for  VOR,  and  V(p,  d)  for  TACAN.  These  empirical 
functions  were  calculated  in  accordance  with  methods  given 
in  NBS  Report  6767  [  1 9 6 1  ]  .  The  V(p,  9)  curves  were 
derived  from  data  mostly  in  the  190  Me/ s  range  and  at 
distances  beyond  the  radio  horizon,  and  were  therefore 
thought  well  applicable  to  the  VOR  frequencies  at  angular 
distances  greater  than  10  milliradians,  whereas  the  V(p,  d) 
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functions  contain  a  frequency  factor,  and  are  based  on  a 
more  extensive  data  sample  extending  up  to  the  frequency 
range  of  TACAN,  and  are  also  more  suitable  for  distances 
within  and  slightly  beyond  the  radio  horizon  than  the  V(p,  9) 
curves.  Thus  the  V(p,  d)  functions  were  used  for  VOR  at 
distances  corresponding  to  0  <  10  milliradians,  and  for 
TACAN  at  all  distances. 

Either  one  of  the  V(p,  0)  or  V(p,  d)  functions  are  used  to 
determine  the  cumulative  distribution  with  time  of  hourly 
median  basic  transmission  loss  values  relative  to  the 
reference  values  calculated  under  (a)  above.  The  results 
for  VOR  and  TACAN  consist  in  expected  distributions  of 
all  hourly  median  basic  transmission  loss  values  during 
the  year  as  functions  of  path  length,  angular  distance, 
terminal  height,  and  carrier  frequency. 

(c)  In  addition  to  the  distributions  of  hourly  medians  represent¬ 
ing  long-term  variations,  short-term  (usually  within-the- 
hour)  distributions  of  the  received  signal  levels  have  to  be 
considered.  Short-term  variations  in  this  particular 
application  are  principally  due  to  two  causes.  One  is  the 
inherent  short-term  fluctuation  of  the  tropospheric  signal 
ascribed  to  the  phase  interference  of  rays  reflected  from 
small  layers  or  scattered  from  refractive  index  discon¬ 
tinuities,  or  to  reflections  from  ground  irregularities.  The 
other  cause  is  the  pattern  of  the  aircraft  antenna:  the 
numerous  small  lobes  and  gain  changes  with  varying  bearings 
are  best  represented  by  a  cumulative  distribution  of  effective 
antenna  gain  with  time  as  the  aircraft  moves  through  space. 
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This  is  actually  a  variation  in  transmission  loss  and  not  in 
basic  transmission  loss  as  it  concerns  the  gain  of  the  receiv¬ 
ing  antenna  in  the  aircraft.  However,  the  actual  calculations 
for  both  VOR  and  TACAN  are  handled  most  efficiently  by 
including  the  gain  of  the  transmitting  antenna  (ground  facility) 
in  the  auxiliary  propagation  curves  developed  for  the  model, 
and  by  treating  the  variable  gain  of  the  receiving  antenna  as 
an  additional  time  variability,  which  has  to  be  combined  with 
the  short-term  fluctuations  due  to  propagation  and  the  long¬ 
term  variability  of  hourly  medians  in  order  to  arrive  at  a 
complete  picture  of  time  variations  which  can  be  used  for 
the  service  and  interference  analysis.  This  will  be  discussed 
in  more  detail  in  the  Appendix,  where  the  mechanics  of  com¬ 
bining  cumulative  distribution  will  be  explained. 

Both  the  calculation  of  service  volumes  limited  by  receiver  noise  only 
and  the  calculation  of  service  volumes  under  various  conditions  of 
interference  require  a  knowledge  of  the  time  distribution  of  transmis¬ 
sion  loss  or  field  strength  at  many  points  in  space;  therefore  time 
variability  calculations  have  to  be  performed  for  these  points  with 
distance  from  the  desired  and  undesired  station,  aircraft  height, 
spacing  of  desired  and  undesired  ground  facilities,  carrier  frequency, 
and  antenna  patterns  as  parameters. 

It  should  be  noted  here  that  the  one -hour  period  taken  as  the 
dividing  line  between  long-term  and  short-term  variations  is  somewhat 
arbitrary.  It  is  convenient  in  view  of  the  available  empirical  time 
variability  functions. 
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3.  INTERFERENCE  BETWEEN  TWO  STATIONS 

As  shown  by  Fig.  1,  both  the  desired  and  the  undesired  signals 
arrive  at  the  aircraft  over  propagation  paths  characterized  by  the 
distances  d  and  d  ,  respectively,  and  by  the  aircraft  height.  The 
distances  are  measured  along  the  great  circle  path.  Both  signals  vary 
with  time,  and  the  distributions  of  signal  levels  are  calculated  in 
accordance  with  the  procedures  outlined  in  the  preceding  section.  The 
next  step  is  the  determination  of  ratios  of  desired  to  undesired  signal 
exceeded  for  given  percentage -of-time  values  when  the  aircraft  is  at  a 
particular  location  in  space.  A  desired  ratio  value  cannot  be  obtained 
directly  for  a  desired  percentage  value  of  time;  thus,  ratio  values  are 
determined  for  a  finite  number  of  points  scattered  over  the  potential 
service  volume  for  9  5%  time  availability,  and  the  location  of  the  three- 
dimensional  9  5%  time  availability  contour  for  a  given  ratio  value  is 
determined  by  interpolation.  This  procedure  has  to  be  repeated  for 
each  set  of  parameters  such  as  spacing  between  co-channel  ground 
facilities,  or  aircraft  altitudes. 

The  ratio  of  de sired-to-unde  sired  signal  can  be  expressed  as  a 
decibel  difference  of  desired  and  undesired  signal  levels,  which  are 
obtained  directly  from  the  calculated  transmission  loss  values  and  the 
system  parameters.  The  distribution  of  the  ratio  will  be  denoted 
D/U(p,  d).  This  notation,  besides  the  separately  specified  aircraft 
height  and  the  spacing  between  the  desired  and  the  undesired  ground 
station,  implies  a  pe rcentage -of-time  value,  p,  as  well  as  a  distance 
parameter,  d,  which  has  to  be  identified  in  each  specific  case.  The 
percentage  value,  p,  is  that  percentage  of  time  during  which  the  given 
value  of  D/U  is  exceeded.  By  virtue  of  the  aircraft  being  in  motion, 
time  variations  also  include  variations  in  space.  As  the  actual  time 
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distribution  of  D/U  may  vary  from  installation  to  installation  due  to 
terrain  characteristics  and  other  factors  not  taken  into  account  in  this 
analysis,  the  time  availability  p  may  be  interpreted  as  an  expression 
of  reliability  for  an  average  installation.  The  concepts  of  prediction 
uncertainty  and  service  probability  in  the  sense  defined  by  Barsis, 
et  al.  [1961]  are  not  used  here  explicitly.  Because  the  calculations 
are  based  on  smooth  earth,  the  standard  deviation  characterizing  a 
prediction  uncertainty  can  be  assumed  to  be  quite  small;  therefore,  the 
expected  value  corresponding  to  a  service  probability  of  0.  5  and  used 
here  is  very  little,  if  any,  different  from  values  applicable  to  a  high 
degree  of  service  probability.  It  should,  however,  be  understood  that 
installations  in  very  rough  terrain  would  show  somewhat  different 
results,  if  analyzed  in  detail.  This  consideration  is  applicable  to  any 
value  of  the  time  availability  p  considered;  the  0.  5  service  probability 
has  no  relation  to  the  percentage  of  time  used,  and  may  approximately 
be  identified  with  the  results  obtained  for  an  average  installation  over 
relatively  smooth  terrain. 

As  an  example,  D/U(95,  d^.)  =  10  db  means  that  for  the  average 
installation  the  ratio  of  the  desired  to  the  undesired  signal  is  equal  to 
or  greater  than  10  db  during  9  5%  of  the  time  at  a  distance  d^  from 
the  interfering  station  for  a  given  aircraft  height  and  spacing  between 
the  desired  and  the  undesired  station. 

In  order  to  solve  for  the  time  availability  of  the  desired-to- 
unde  sired  ratios  at  any  point  in  space,  time  distributions  of  the  indi¬ 
vidual  power  or  signal  levels  (expressed  in  decibels)  have  to  be  com¬ 
bined  in  a  similar  way  as  was  done  to  obtain  the  total  time  variability 
of  transmission  loss.  Thus  a  time  distribution  of  the  decibel  differ¬ 
ences  (representing  the  power  ratios)  can  be  obtained,  from  which  the 
time  availability  of  the  specified  ratio  value  can  be  determined  for 


each  point  in  space.  The  method  of  combining  cumulative  distributions 
will  be  illustrated  in  the  Appendix,  as  mentioned  before. 


4.  SYSTEM  PARAMETERS 

System  parameters  are  based  on  a  standard  VORTAC  ground 
installation,  and  passenger-jet  aircraft  antennas.  A  52  ft  diameter 
counterpoise,  1  2  ft  above  ground,  is  used  at  the  ground  station  for 
both  the  VOR  and  the  TACAN  antenna.  Other  parameters  are  as 
follows: 

TABLE  I 

System  Parameters 


Carrier  Frequency  Used  in 
the  Calculations 

Transmitter  Power 

Transmission  Line  Losses 

Effective  Radiated  Power 

Polarization 

Ground  Antenna  Type 

Maximum  Antenna  Gain 
Relative  to  Isotropic  Antenna 

Horizontal  Pattern 

Vertical  Pattern 

Antenna  Height  above 
Counte  rpoise 

Antenna  Height  above  Ground 

Aircraft  Antenna  Type 

Aircraft  Antenna  Gain 

Required  Receiver  Input 
Class  I 
Class  II 


VOR 

TACAN 

113  Mc/s 

1150  Me  / s 

23  dbw 

- 

6  db 

3  db 

- 

39,  43,  &  48  dbw 

Horizontal 

Ve  rtical 

F  our-loop 

Alford  Array 

Center  Array 

2.  15  db 

8.  15  db 

Appr.  Circular 

Appr.  Circular 

Appr.  like  Dipole 

See  Fig.  3 

4  ft 

18  ft 

1  6  ft 

30  ft 

E  -Cavity 

Annular  Slot 

See  Fig.  2 

See  Fig.  4 

5pV  into  50  ohms 

-109  dbw 

OpV  into  50  ohms 

- 1 03  dbw 

-11- 


For  the  VOR  system,  the  6  db  line  loss  figure  includes  trans¬ 
mitting  and  receiving  transmission  line  and  associated  losses.  The 
gain  of  the  aircraft  antenna  was  determined  as  a  function  of  azimuth 
and  vertical  angle  from  a  modeling  study  based  on  an  E -cavity  type 
VOR  antenna  in  the  vertical  stabilizer  of  a  pas senge r-type  jet  aircraft 
[Conviar,  1959]  .  Fig.  2  shows  the  distribution  of  antenna  gains  at 
various  vertical  angles,  together  with  the  measurement  points  from 
which  the  principal  distribution  was  derived. 

The  several  values  of  effective  radiated  power  assumed  for  the 
TACAN  model  include  the  line  losses  in  the  ground  station  and  the 
transmitter  power;  therefore,  the  line  loss  figures  given  in  Table  I 
are  the  losses  in  the  aircraft  receiving  system.  The  ground  antenna 
was  taken  to  be  the  typical  center  array  shown  by  Casabona  [19  56]  , 
with  its  free -space  vertical  radiation  pattern  shown  on  Fig.  3.  The 
gain  of  the  aircraft  antenna  again  was  obtained  from  modeling  studies 
based  on  an  annular  slot  type  DME  antenna  mounted  on  the  bottom 
fuselage  center  line  forward  of  the  landing  gear  of  a  passenger  jet 
aircraft.  Pertinent  distributions  of  antenna  gain  values  for  several 
vertical  angles  are  shown  on  Fig.  4,  together  with  the  results  of 
available  measurements  [Commercial  Jetstar,  19  59]  . 

5.  COMPUTER  PROGRAMS 

Calculations  necessary  to  produce  all  desired  distributions  of 
transmission  loss  values  or  power  levels  as  well  as  the  processes 
involving  the  combination  of  various  distributions  are  so  numerous 
and  complex  that  the  use  of  a  large  electronic  computer  is  mandatory. 
Computer  programs  were  organized  in  the  following  manner. 


-12- 


(a)  Existing  programs  were  modified  and  new  programs  were 
developed  to  calculate  reference  basic  transmission  loss 
values  including  the  effect  of  the  transmitting  antenna  as  a 
function  of  distance  and  aircraft  height  for  both,  VOR  and 
TACAN  separately.  The  resulting  propagation  curves  were 
reduced  to  tabular  form,  and  made  available  to  the  com¬ 
puter  as  "data"  for  subsequent  procedure. 

(b)  Subroutines  were  developed  for  the  combination  of  various 
cumulative  time  distributions.  These  served  either  to 
combine  the  various  variability  distributions,  or,  with 
changed  sign,  to  obtain  the  final  result  of  the  computer 
programs,  namely  the  distribution  of  de sired-to-unde sired 
signal  ratios . 

(c)  Short-term  and  long-term  time  distribution  functions  were 
made  available  to  the  computer  in  tabular  form  as  "data" 
separately  for  VOR  and  TACAN.  Distributions  of  receiving 
antenna  gain  values  were  handled  in  a  similar  way. 

(d)  The  main  computer  program  used  all  the  "data"  assembled 
by  the  subroutines  described,  and  delivered  tables  of  inter¬ 
ference  ratios  for  various  values  of  time  availability  as  a 
function  of  the  aircraft  position  and  various  distances  between 
the  desired  and  the  undesired  facilities.  From  these  tables, 
auxiliary  graphs  were  plotted,  and  the  9  5%  time  availability 
values  determined  which  are  the  results  of  the  study,  and 
will  be  shown  in  graphical  form  in  the  next  section. 
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The  final  program  step,  after  check-out  and  after  assembly  of 
all  subroutines  and  all  "data,  "  consumed  about  15  hours  of  computer 
time.  It  is  not  possible  to  estimate  the  number  of  man  hours  required 
to  accomplish  the  same  object  by  the  use  of  slide  rule  and  desk 
calculator . 


6.  RESULTS  OF  THE  STUDY 

The  results  are  in  the  form  of  prediction  curves.  Three  types 
of  curves  are  presented,  which  involve  service  limitations  due  to: 

(1)  receiver  sensitivity  and  transmitter  power  output  only,  i.e.,  the 
service  is  limited  by  noise;  (Z)  interference  from  one  co-channel 
station;  and  (3)  proximity  to  an  undesired  station  which  may  be  operat¬ 
ing  on  an  adjacent  channel. 

a.  Service  Limitations  Due  to  Noise 

Curves  of  the  first  type  are  shown  as  solid  lines  on  Figs.  5  and 
6  labeled  "VORTAC  (or  TACAN)  Service  Volume  without  Interference" 
and  as  dashed  lines  on  Figs.  7  through  12  labeled  "VOR  (or  TACAN) 
Service  Volumes  with  Interference  from  One  Co-Channel  Station.  " 

In  the  volume  defined  by  the  revolution  of  the  appropriate  curve  about 
the  ordinate  axis  (the  desired  station  is  near  the  origin)  service  as 
defined  in  Section  3  is  predicted  to  be  available  during  9  5%  of  the  time. 
This  may  be  called  a  conservative  estimate,  as  interference  conditions 
are  based  on  the  values  obtained  on  the  direct  path  between  the  desired 
and  the  undesired  station;  for  other  points  on  the  surface  of  revolution 
there  is  obviously  less  interfering  signal.  The  appropriate  curve  is 
determined  by:  (1)  system  type  (VOR  or  TACAN),  (2)  receiver  sensi¬ 
tivity  (Class  I  or  Class  II),  and  (3)  ground  station  power  (VOR,  23  dbw 
transmitter  power;  TACAN,  48,  43,  or  39  dbw  effective  radiated  power). 
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The  system  parameters  have  been  described  in  Section  4  above; 
however,  the  ground  station  power  output  for  the  "standard"  TACAN 
curves  on  Fig.  5  was  taken  as  the  48  dbw  E.R.P.  (effective  radiated 
power)  value.  The  effect  of  changes  in  the  values  of  the  E.R.P.  are 
illustrated  on  Fig.  6,  which  shows  Class  I  and  Class  II  TACAN  service 
volumes  for  the  additional  39  and  43  dbw  E.R.P.  values,  constituting 
noise -limited  service  volumes. 

The  Class  I  and  Class  II  noise -limit  curves  for  VOR  and  standard 
TACAN  ( 48  dbw  E.R.P.  )  are  repeated  on  Figs .  7  through  9  and  1 0 
through  12,  respectively  in  order  to  show  where  noise  is  more 
restrictive  than  co-channel  interference. 

All  figures  discussed  so  far  show  also  the  distance  to  the  radio 
horizon  as  a  function  of  aircraft  altitude,  based  on  the  antenna  heights 
over  a  smooth  earth  and  a  standard  atmosphere.  Due  to  the  higher 
elevation  of  the  TACAN  antenna,  its  radio  horizon  extends  about  two 
miles  further  for  each  aircraft  altitude  than  the  VOR  horizon. 

b.  Service  Limitation  Due  to  Co-Channel  Interference 

The  second  type  of  curves,  shown  on  Figs.  7,  8,  and  9  for  VOR 
and  10,  11,  and  12  for  TACAN,  illustrates  the  effect  of  co-channel 
interference  on  the  service  volume  when  the  aircraft  is  located  above 
the  great  circle  path  between  the  desired  and  the  undesired  station  at 
a  distance  d  from  the  desired  station.  The  geometry  involved  is 
shown  by  a  small  diagram  on  each  curve  sheet.  Station  separations, 

S,  ranging  from  70  to  695  nautical  miles  were  considered  along  with 
aircraft  altitudes  ranging  from  1,  000  to  100,  000  feet.  Each  curve 
sheet  is  applicable  to  a  different  signal-to-interference  ratio, 

D/U(95,  d).  For  example,  D/U(95,  d)=  14  db  means  that  the  desired 
signal  is  at  least  14  db  greater  for  9  5%  of  the  time  along  the  solid 
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curve  which  forms  the  boundary  of  service  volume  closest  to  the  inter¬ 
fering  station.  On  these  figures  the  limitation  imposed  by  ground 
station  power  output  and  receiver  sensitivity  is  described  only  by  the 
dashed  lines  labeled  Class  I  or  Class  II. 

The  volume  defined  by  rotating  the  appropriate  curve  about  the 
ordinate  axis  represents  a  volume  in  which  service  reliability  in  the 
sense  defined  in  Section  3  is  9  5%  or  greater  (by  virtue  of  each  curve 
representing  the  smallest  station  separation).  Similarly,  if  service 
is  limited  by  interference  from  several  co-channel  stations,  the 
volume  defined  by  rotating  the  most  restrictive  curve  (which  is  the 
curve  appropriate  to  the  closest  interfering  station)  about  the  ordinate 
axis  represents  a  volume  in  which  the  service  reliability  is  generally 
95%  or  greater,  when  an  average  installation  (in  the  sense  described 
above)  is  considered. 

Figures  7,  8,  and  9  describe  VOR  service  volumes  with  one 
co-channel  interfering  station  for  signal-to -inte rfe rence  ratios  of 
D/ U  (95,  d)  =  14,  20,  and  26,  re  spectively . 

Figures  10,  11,  and  12  describe  TACAN  service  volumes  with 
one  co-channel  interfering  station  for  signal-to -inte  rfe  rence  ratios 
of  D/ U  (9  5,  d)  =  8,  14,  and  20  db,  re  spectively. 

c.  Signal-to-Interference  Ratio  Curves 

The  third  type  of  curves  shows  signal -to -inte  rference  ratios  in 
the  neighborhood  of  an  undesired  station.  The  distance,  d^,  from 
the  undesired  station  is  used  to  locate  the  aircraft  in  the  same  manner 
as  d  was  used  for  the  service  volume  curves  in  the  preceding  figures. 
This  geometry  is  shown  by  the  small  diagram  on  each  curve  sheet. 
Signal-to-interference  ratio  values  of  D/U(95,  d)  are  used  as  the 
ordinate,  with  emphasis  on  negative  values.  Curves  applicable  to 
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station  separations  ranging  from  25  to  175  nautical  miles  are  presented 
on  some  of  the  graphs.  There  is  a  separate  figure  for  each  system- 
altitude  combination  with  the  altitudes  ranging  from  1,  000  to  100,  000 
feet.  As  an  example,  the  intersection  of  the  D/U(95,  d)  =  -20  db 
ordinate  on  Fig.  13  with  the  curve  labeled  S  =  70  means  that  for  95% 
of  the  time  the  interfering  signal  exceeds  the  desired  signal  by  at 
least  20  db  at  a  distance  of  27  miles  from  the  undesired  station  and 
for  the  1,  000  ft  aircraft  height  represented  by  Fig.  13.  If,  in  an 
application  to  an  adjacent  channel  interference  problem,  the  undesired 
signal  just  causes  trouble  when  it  is  greater  than  the  desired  signal 
by  20  db,  service  is  available  95%  of  the  time  where  the  D/U(95,  d) 
curve  for  the  appropriate  station  separation  does  not  become  more 
negative  than  -20  db. 

The  curve  sheets  developed  for  signal  ratios  near  an  undesired 
station  are  arranged  as  shown  on  Table  II: 


TABLE  II 

Index  of  Curve  Sheets 


Figure 

Numbe  r 

Aircraft  Altitude 

VOR 

TACAN 

In  Feet 

13 

26 

1,  000 

14 

27 

5,  000 

15 

28 

10, 000 

16 

29 

15,  000 

17 

30 

20,  000 

18 

31 

30, 000 

19 

32 

40,  000 

20 

33 

50,  000 

21 

34 

60,  000 

22 

35 

70,  000 

23 

36 

80,  000 

24 

37 

90,  000 

25 

38 

100, 000 
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The  curves  on  Figs.  13-38  do  not  show  the  effect  of  interference 
beyond  the  undesired  station  or  at  locations  not  on  the  great  circle  path 
connecting  the  stations;  i.e.  the  distance  d^.  shown  on  the  abscissa 
scales  locates  the  aircraft  on  the  great  circle  path  between  the  desired 
and  the  undesired  station.  However,  a  method  of  approximating  the 
locus  of  a  constant  interference  ratio  D/U(95,  d^)  as  a  circle  about 
the  undesired  station  has  been  developed.  For  a  given  aircraft  altitude 
this  circle  is  centered  on  the  extension  of  the  line  connecting  the  ground 
stations  concerned,  but  on  the  "far"  side  of  the  undesired  station.  The 
pertinent  geometry  is  shown  in  a  top  and  a  side  view  on  Fig.  39.  Gen¬ 
erally  service  may  be  regarded  as  being  unsatisfactory  within  this 
circle,  even  though  some  locations  having  satisfactory  service  may 
exist  above  the  undesired  station  due  to  the  vertical  pattern  of  the 
ground  antenna. 

Two  basic  assumptions  have  to  be  made  as  follows: 

(a)  The  geometry  represented  by  Fig.  39  may  be  considered 
plane  geometry;  thus  the  earth  is  assumed  to  be  flat  and 
the  projection  of  the  slant  ranges  r  and  r  onto  the 
horizontal  plane  is  approximately  equal  to  the  actual 
ranges  r^  and  r^. 

(b)  The  interference  ratio  D/U(95,  d  )  is  proportional  to 
the  logarithm  of  the  ratio  of  the  ranges  rjj/  ry*  In  the 
terminology  used  in  the  appendix,  this  may  be  written 
for  any  percentage  value  p: 

D/U(p,  du)  =  L(p,  dy)  *  [-L(p,  dD]  =  Mlogl^/^)  (1) 


where  M  is  a  constant. 
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Assumption  (a)  is  reasonable  if  A  (aircraft  altitude/ 1 6 00), 

where  the  distances  d^  and  d  are  in  nautical  miles  and  the  aircraft 
altitude  is  in  feet.  Hence,  the  problem  of  finding  the  locus  of  r^/ r^ 

=  constant  may  be  solved  using  plane  geometry.  This  locus  is  found 
to  be  a  circle,  described  by  the  parameters  C  and  R  defined  on 
Fig.  39,  and  given  by  the  following  equations  which  contain  an  addi¬ 
tional  auxiliary  parameter  B: 


B 


d  /d  =  (S  -  d  )/ d 
D  U  U  U 


(2) 


S(BZ  +  1) 

2(B2  -  1) 


R 


(3) 

(4) 


As  an  example,  use  the  separation  S  =  1  50  miles  for  the  VOR 
graph  of  Fig.  1  6,  representing  an  aircraft  height  of  15,  000  ft.  For  a 
desired-to-undesired  signal  ratio  D/U(95,  d  )  =  -40  db  which  may 
very  well  constitute  the  adjacent-channel  interference  threshold  in  a 
practical  case,  the  distance  d^.  =  20  miles  from  the  undesired  station 
(with  d^  =  130  miles)  is  read  from  Fig.  16.  The  results  are: 


B 


130 

20 


6 .  5 


C  =  75  (43/41)  =  79  miles 
R  =  97  5/41  =  24  miles 
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Thus,  inside  a  circle  of  24  miles  radius  centered  at  a  point  4  miles 
beyond  the  undesired  station,  VOR  service  at  15,  000  ft  aircraft  alti¬ 
tude  would  be  expected  to  be  sub-standard. 

It  should  not  be  overlooked  that  this  method  is  only  approximate. 
The  second  assumption,  (b),  would  be  violated  by  lobing  in  the  trans¬ 
mission  loss  versus  distance  curve,  which  may  occur  at  any  constant 
altitude  due  to  ground  reflections  or  due  to  the  ground  antenna  pattern. 
The  assumption  (b)  would  also  be  violated  by  a  change  in  the  slope  of 
the  transmission  loss  versus  distance  curve  which,  as  an  example, 
occurs  in  the  vicinity  of  the  radio  horizon.  However,  the  use  of  the 
large st  d  read  for  a  particular  station  spacing  and  a  particular 
signal  ratio  from  Figs.  13-38  avoids  the  ambiguity  due  to  lobing. 
Furthermore,  in  most  applications  service  is  limited  by  noise  rather 
than  by  interference  if  ranges  beyond  the  radio  horizon  of  the  desired 
station  are  considered. 

7.  COMPARISON  OF  THE  RESULTS  WITH  PREVIOUS  STUDIES 

Staras,  Rice,  and  Herbstreit  [1951]  have  previously  predicted 
VOR  service  volumes.  The  curves  they  developed  are  not  in  overall 
agreement  with  curves  presented  here.  Although  the  antennas,  both 
ground  and  air,  used  in  the  two  predictions  were  not  identical,  they 
were  so  similar  that  only  minor  differences  between  prediction  curves 
would  be  likely  because  of  different  antennas.  However,  since  19  51 
more  refined  methods  for  calculating  transmission  loss  beyond  the 
radio  horizon  and  for  estimating  time  variability  have  been  developed. 
Thus,  major  differences  in  prediction  curves  are  most  likely  caused 
by  different  variability  estimates,  and  different  predicted  median  fields 
beyond  the  radio  horizon.  The  power  output  for  VOR  used  in  the 
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earlier  study  was  3  db  lower  than  that  used  here,  but  only  the  curves 
for  service  without  interference  would  be  affected  by  this  difference. 

The  TACAN  curves  developed  by  Decker  [1956]  also  differ  from 
the  curves  presented  here.  The  same  ground  station  antenna  was  used 
in  both  studies,  but  in  the  earlier  study  it  was  placed  100  feet  above 
the  ground,  whereas  in  this  study  a  height  of  30  feet  was  used.  Although 
the  method  used  to  account  for  ground  reflections  was  almost  identical, 
the  actual  numbers  involved  were  different,  because  of  the  different 
ground  antenna  height  and  different  estimates  of  the  effective  ground 
reflection  coefficient.  This  may  account  for  minor  differences  in  the 
prediction  curves,  but  different  estimates  of  long-term  variability, 

V(p,  d),  were  found  to  be  more  significant.  Decker's  probability  of 
service  curves  applicable  to  interference  limited  service  also  included 
limitations  due  to  noise  (transmitter  power  and  receiver  sensitivity), 
whereas  in  the  present  study  interference  and  noise  limitations  were 
treated  separately. 
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SKETCH  SHOWING  RELATIVE  POSITIONS  OF  AIRCRAFT  AND 
OF  NAVIGATIONAL  FACILITIES  OVER  A  SMOOTH  SPHERICAL  EARTH 
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CUMULATIVE  DISTRIBUTION  OF  AIRCRAFT  ANTENNA  GAIN 
VOR,  STABILIZER  MOUNTED  E-CAVITY,  HORIZONTAL  POLARIZATION,  II3MC/S 
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TYPICAL  URN-3  ANTENNA  PATTERN 

Free  Space  Radiation  in  the  Vertical  Plane 


I 


O  CO  -'O  ^  rv] 

O  O  O  O 


T-* 


TO 


to 

W 

W 

2 

a 

w 

Q 

2 


u 

0 


o 

>- 

W 

2 

o 

2 

< 

2 

0 
I — l 

H 

C 

> 

W 

2 

W 


ro 

3 

oc 

•  i-l 
2 


8  jo 


§  MIVD  30VI10A  3AIXV33TT 


B-  33540 


CUMULATIVE  DISTRIBUTION  OF  AIRCRAFT  ANTENNA  GAIN 

TACAN  ,  BELLY  MOUNTED  ANNULAR  SLOT,  VERTICAL  POLARIZATION,  1150  Mc/S 


PERCENT  OF  OCCURRENCES  IN  WHICH  GAIN  EXCEEDS  ORDINATE  VALUE 


VORTAC  SERVICE  VOLUME  WITHOUT  INTERFERENCE 
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TACAN  SERVICE  VOLUME  WITHOUT  INTERFERENCE 
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APPENDIX 


1.  VOR  PROPAGATION  MODEL 

As  mentioned  briefly  in  Section  2  of  the  main  body  of  this  report, 
the  smooth,  spherical  earth  model  was  used  to  calculate  transmission 
loss  for  the  VOR  system  at  113  Mc/s.  A  linear-gradient  atmosphere 
was  assumed  for  the  initial  calculations,  so  that  for  a  first  approxi¬ 
mation  radio  rays  could  be  considered  to  be  straight  lines  above  an 
earth  having  an  effective  radius  k  times  its  actual  value.  Where 
appropriate,  ray  tracing  methods  were  used  for  final  calculations.  In 
conformance  with  usual  practice,  the  effective  earth  radius  factor,  k, 
was  assumed  to  be  4/3  for  the  "standard"  atmosphere.  Calculation 
methods  are  based  on  material  contained  in  NBS  Report  67  67  [1961]  . 

It  was  found  convenient  to  include  only  the  gain  of  the  transmitting 
antenna  in  the  definition  of  transmission  loss  as  used  in  these  calcula¬ 
tions.  Usually,  the  term  transmission  loss  includes  the  antenna  gains 
at  both  path  terminals,  whereas  the  term  "basic  transmission  loss" 
excludes  all  antenna  gains  [Norton,  19  53  and  19  59]  .  In  this  study,  an 

hourly  median  transmission  loss,  L  ,  is  defined  to  include  the  free- 

m 

space  gain  of  the  transmitting  antenna  as  a  function  of  the  elevation 
angle  above  ground  as  well  as  the  effect  of  reflections  from  the  counter¬ 
poise  and  from  the  ground.  The  transmitting  antenna  at  the  ground 
VOR  station  is  a  standard  horizontally  polarized  Alford  loop  array  4  ft 
above  the  center  of  a  52  ft  diameter  counterpoise,  which  in  turn  is  1  2  ft 
above  ground  level. 

Within  the  radio  horizon,  values  of  the  hourly  median  transmis¬ 
sion  loss,  L  ,  were  calculated  using  geometrical  optics  methods 
m 

[Kirby,  Herbstreit,  and  Norton,  1952]  .  These  methods  take  into 


1-2 


account  the  interference  between  the  direct  and  the  ground-reflected 
ray.  Fig.  1-1  shows  the  geometry  for  this  ray  interference  problem, 
and  defines  many  of  the  symbols  used  in  the  analysis.  Using  this  geo¬ 
metry,  two  curves  of  L  versus  distance  were  calculated  for  each 

m 

assumed  aircraft  altitude.  One  curve  applies  to  reflections  from  the 
counterpoise,  where  applicable,  and  the  other  to  reflections  from  the 
ground.  Both  curves  were  blended  in  the  region  where  the  assumed 
specular  reflection  point  is  close  to  the  edge  of  the  counterpoise.  In 
order  to  facilitate  this  blending,  the  curves  were  calculated  so  that 
they  extended  beyond  their  region  of  validity. 

For  reflections  from  the  counterpoise,  calculated  for  elevation 
angles  greater  than  5  degrees,  the  effective  earth  radius  is  increased 
by  the  height  of  the  counterpoise  above  ground  so  that  a  sphere  through 
the  level  of  the  counterpoise  becomes  the  reflecting  surface.  For  the 
calculations,  all  heights  and  distances  are  expressed  in  statute  miles, 
and  all  angles  in  radians.  The  following  constants  and  relationships 
were  used,  with  most  of  the  symbols  defined  on  Fig.  1-1. 

a  .  .  . .  effective  earth  radius  =4/3  times  the  actual 

radius  of  the  earth  (a  =  5280  statute  miles). 


a 

c 


h 


1 


h 


2 


effective  earth  radius  for  reflections  from, 
the  counterpoise  (a  =  a  +  1 2  ft,  but  is  also 
expressed  in  statute  miles). 

height  of  the  transmitting  antenna  above  the 
reflecting  surface  (h^  =  4  ft  for  reflections 
from  the  counterpoise,  and  1  6  ft  for  reflec¬ 
tions  from  the  ground). 

aircraft  altitude  above  the  reflecting  surface 
( h 2  is  the  aircraft  altitude  above  ground  for 
reflections  from  the  ground,  and  the  altitude 
minus  1  2  ft  for  reflections  from  the  counter¬ 
poise). 
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|R| .  magnitude  of  the  reflection  coefficient,  assumed 

to  be  0.95  for  counterpoise  reflections  and  0.9 
for  ground  reflections. 

c . .  phase  angle  of  the  reflection  coefficient  relative 

to  tt  ,  and  assumed  to  be  zero  for  both  cases. 


D. .  divergence  factor,  assumed  to  be  unity  for 

reflections  from  the  counterpoise. 


d .  distance  along  the  great  circle  path  to  a  point 

on  the  reflecting  surface  below  the  aircraft. 


r 

o 


length  of  the  direct  ray  path. 


r1+r2.... 


length  of  the  reflected  ray  path. 


d 


1 


distance  from  transmitting  antenna  to  reflection 
point. 


d^  .......  distance  from  aircraft  to  reflection  point. 

iJj .  grazing  angle  at  reflection  point. 


The  first  step  was  the  determination  of  the  grazing  angle  l|j.  The  con¬ 
stant  a  shown  in  the  equations  below  refers  to  reflections  from  the 
ground  (i \j  <  15°);  for  reflections  from  the  counterpoise  (iji  >  5°),  the 
effective  earth  radius  a^  was  used  instead  of  a,  and  the  heights  h^ 
and  h^  were  changed  accordingly.  Initial  values  of  the  linear  dis¬ 
tance  d  are  estimates,  which  were  checked  later  by  ray-tracing 
methods,  and  changed  if  necessary: 


d  =  d  +  d2  (I- 1 ) 

h  1  =  h  -  d  Z/2a  (I-2a) 

h  1  =  h  -  d  2/2a 


(I-2b) 
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tan  +  =  hj'/dj  =  h2'/d2  (1-3) 

Equations  (I-Za),  (I-Zb),  and  (1-3)  are  approximate  relations.  In 
order  to  utilize  the  electronic  computer  in  solving  the  above  expressions 
for  the  grazing  angle  i|j,  an  iteration  method  was  used,  assuming  for 
each  d  different  values  of  d  until  the  equations  were  satisfied  within 
the  desired  limits  of  accuracy,  which  was  ±10  miles  in  this  case. 
Final  values  for  ijj  for  each  assumed  total  distance  d  were  then  used 
to  calculate  various  auxiliary  angles  and  quantities  to  be  used  in  the 
ray-tracing  problem  and  in  the  final  calculation  of  the  transmission 

loss  L  .  For  reflections  from  the  counterpoise  h  '  may  be  assumed 

m  1 

equal  to  h ^ ,  d?  to  equal  d,  and  iJj  can  be  determined  directly  with¬ 
out  iteration. 

Referring  to  Fig.  1-1,  the  ray  path  portion  r  for  reflections 
from  the  counterpoise  can  be  very  closely  approximated  by: 

ri  ~  hl// Sin  U"4) 

as  the  distance  to  the  reflection  point  is  extremely  small.  By  the  same 
reasoning, 

cos  ^)/ac  (1-5) 

Other  needed  quantities  are  calculated  as  follows: 


sin  4>  =  (a  cos  dj)/(a  +  h  ) 

Z  c  c  Z 


d-6) 


e2  =  tt/2-4j-4>z 


(1-7) 
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e  =  e  +  e_ 

o  12 

r2  =  (ac  Sin  02^^sin  $2) 

r  =  Mh  -h)Z+4(a  +  h  )  (a  +  h  )  sin  (9/2)  (I- 10) 

o  \J  1  2  clc2  o 

The  total  angular  distance  0^  is  now  used  to  obtain  the  corrected 
great-circle  distance  d  from 

d  =  a  0  ,  (I- 11) 

o 

which  may  differ  slightly  from  the  first  assumption  given  by  (1-1). 

For  the  reflections  from  the  ground,  the  approximations  (1-4) 
and  (1-5)  can  no  longer  be  used,  because  the  reflection  point  may  not 
be  close  enough  to  the  transmitting  antenna.  Referring  again  to  Fig. 

I  —  1 ,  the  following  relations  are  established,  again  after  first  obtaining 
the  grazing  angle  i]j  : 


(1-8) 

d-9) 


sin  <|>  = 

P 

o 

o 

03 

PT 

+ 

l — 1 

d-12) 

sin  c|>  = 

(a  cos  \\>)/ (a  +  h_,) 

d-13) 

0 

1 

=  tt/z-4j-<4>i 

d-14) 

h 

=  tt  /  Z  -  -  *)>2 

(1-15) 

rl  = 

(a  sin  0  )/( sin  <)>  ) 

(1-16) 
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r7  =  (a  sin  9  )/(sin  <f>2) 


d-17) 


9  =  0  +  0 

o  1  Z 


d-8) 


r 

o 


+  4(a  +  h  )  (a  +  h  )  sin  (9  / Z) 
1  Z  o 


(1-18) 


Again  the  corrected  great  circle  distance  d  is  obtained  using  ( I—  1 1 ) . 
For  this  case  of  reflections  from  the  ground,  a  divergence  factor  D 
is  now  calculated  using: 


D  = 


a(ri  +  r^)  cosi)j 
(a  +  h^  )(a  +  h^) 


sin  ijj 


sin  9  (sin  9,  cos  c(>_  +  sin  9_,  cos  (b,  ) 
o  1  Z  Z  1 


(1-19) 


Before  calculating  the  transmission  loss  values  as  a  function  of 
distance  for  reflections  from  the  counterpoise  (t|»  >  5°)  as  well  as  for 
reflections  from  the  ground  (i|j  <  15°),  it  is  necessary  to  determine 
the  angles  y  and  y  ^  for  which  the  voltage  gain  of  the  ground 
antenna  is  used.  From  Fig.  1-1, 


sin  (3  =  [r  sin  (tt  -  Zi|>)] /r 


d-20) 


Y 


o  1 


(3  -  ijj  -  0} 


d-21) 


Y 


rl 


-(^  +  9X) 


(I-2Z) 


1-7 


The  voltage  gains  g^  and  g  of  the  ground  antenna  relative 
to  that  of  an  isotropic  radiator  are  given  by  1 .  28  times  the  cosine  of 
the  elevation  angle  of  the  direct  and  the  ground- reflected  rays, 
respectively: 

g  =  1 . 28  cos  y  (I- 2 3a) 

o  I  o  1 

g  =  1.28  cos  y  (I-23b) 

o2  o2 

All  quantities  necessary  in  the  calculation  of  L  have  now  been 

m 

determined  (r  ,  r,  r  ,  |  R I ,  c,  g  1 ,  g  _,  D,  and  d),  except  the 

o  1  2  ol  o2  r 

relative  phase  angle  between  the  direct  ray  r  and  the  ground- 
reflected  ray,  r ^ .  This  phase  angle  A,  in  radians,  can  be  expressed 
by: 


A 


33.73  f 


Me  1 


( r  +  r 


-  ro>- 


(1-24) 


where  f  =  113  Mc/s  is  the  carrier  frequency  of  the  signal,  taken 
at  the  center  of  the  VOR  band. 

Finally,  the  formula  for  L  can  be  written  as  follows,  where 

m 

the  actual  lengths  of  the  ray  paths  r  ,  r^,  and  r_>  are  used  instead  of 
the  great  circle  distance  d.  The  great  circle  distance,  however,  is 
used  as  an  abscissa  parameter  in  plotting  L  ,  and  for  subsequent 
applications . 


L  =  36.  58  +  20  log  f  +  20  log  r 
m  1 0  Me  1 0  ( 


1 0  log  ^ 
s10 


«ol  + 


r  D  R  g  1 
o  1  1  °rl 


b  +r2 


2r  D  |r|  g  ,g  ,  cosA 
o  1  1  5ol  &r  1 


rl  +r2 


(1-25) 
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For  large  distances  and  low  angles  the  linear  gradient  atmos¬ 
phere  predicts  too  much  bending  [Bean  and  Thayer,  1959]  .  To  correct 
for  this,  a  ray  leaving  the  ground  station  antenna  at  the  same  angle, 
y  as  the  direct  ray  calculated  using  a  linear  gradient  atmosphere 
was  traced  through  an  exponential  atmosphere  until  it  reached  the  air¬ 
craft  altitude.  The  great  circle  distance  below  the  ray  was  then  used 
in  the  final  calculations  instead  of  the  distance  d  obtained  by  (I- 11). 

The  angle,  y  at  which  the  direct  ray  arrived  at  the  aircraft  was 
also  determined  by  this  ray  tracing  procedure.  The  exponential 
atmosphere  used  corresponds  to  the  reference  atmosphere  given  by 
Rice,  Longley,  and  Norton  [19  59]  with  a  surface  refractivity  of 

N  =  301. 
s 

Calculations  of  transmission  loss  beyond  the  radio  horizon  were 
made  by  combining  smooth-earth  diffraction  with  tropospheric  scatter. 
The  short-cut  method  of  Vogler  [1961]  was  used  for  diffraction  compu¬ 
tations,  as  it  specifically  applies  to  horizontal  polarization.  Scatter 
computations  and  methods  of  combining  the  two  mechanisms  follow  the 
work  by  Rice,  et  al.  [1959]  . 

Calculations  for  the  beyond-the -horizon  case  are  based  on  the 
same  exponential  atmosphere  used  for  within-the -horizon  computations. 
For  the  beyond-the -horizon  case  the  effect  of  the  finite  counterpoise 
can  be  neglected,  as  the  grazing  angles  involved  are  very  low  or  zero; 
thus  the  transmitting  antenna  was  treated  as  a  single  loop  1  6  ft  above 
ground.  The  electrical  constants  of  the  ground  assumed  for  clacula- 
tions  of  the  diffracted  wave  were  a  dielectric  constant,  e  =  15,  and 
a  conductivity  value,  cr  =  5  millimhos  per  meter.  At  the  carrier 
frequency  used,  however,  the  effect  of  ground  constants  is  small. 
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Finally,  from  a  large  number  of  computations  both  within  and 
beyond  the  radio  horizon,  a  set  of  curves  was  drawn  for  discrete  air¬ 
craft  altitudes  showing  the  expected  median  transmission  loss  as  a 
function  of  the  great  circle  distance,  d,  and  assuming  the  antenna 
gain  of  the  ground  station  to  correspond  to  a  half-wave  dipole.  The 
construction  of  these  curves  involved  a  certain  amount  of  blending 
from  one  propagation  mechanism  to  another. 

Long-term  variability,  short-term  fading,  and  the  effective  gain 
variations  of  the  airborne  antenna  were  treated  statistically.  Long¬ 
term  variability  is  described  by  the  time  availability  function  V(p,  d) 
for  0  <  10  milliradians  and  V(p,  9)  for  9  >  10  milliradians .  The 
angular  distance,  9,  is  defined  to  be  the  angle  between  horizon  rays 
from  the  transmitter  and  receiver  [Norton,  Rice,  and  Vogler,  1955], 
and  is  defined  on  Fig.  1  of  the  main  body  of  this  report.  For  within- 
the-horizon  cases,  0  is  negative,  and  not  explicitly  used.  Equations 
and  curves  from  NBS  Report  6767  [1961]  along  with  the  curve  by  Rice, 
et  al.  [1959]  giving  effective  heights  for  high  antennas  were  used  to 
calculate  the  V(p,  d)  and  V(p,  0)  functions  applicable  to  all  hours  of 
the  year. 

Short-term  fading  was  described  by  a  cumulative  distribution 

function  V  (p,  9);  it  includes  the  summing  of  a  constant  vector  and  a 
F 

Rayleigh-distributed  vector.  An  auxiliary  distribution  function 
V'  (p,  9)  was  developed  using,  (1)  fading  range  data,  R(0.  1)  -  R(0.9), 
for  100  Mc/s  from  a  paper  by  Janes  [1955]  to  express  fading  range  as 
a  function  of  9,  and  (2)  Fig.  6  and  Table  I  from  Norton,  Vogler, 
Mansfield,  and  Short  [19  5  5]  to  specify  a  cumulative  distribution, 

V'  (p,  0),  for  the  fading  range  given  by  a  particular  value  of  9.  For 
any  percentage  level  R(\/l00)  is  the  magnitude  level  of  a  cumula¬ 
tive  distribution  representing  the  sum  of  a  constant  vector  and  a 
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Rayleigh-distributed  vector  expressed  in  decibels  relative  to  the  con¬ 
stant  component.  Fading  range  as  a  function  of  the  angular  distance, 
9,  in  degrees,  was  approximated  by 


R(0.  1)  -  R(0.  9) 

-  0  for  9  <  0, 

(I-26a) 

R(  0.  1)  -  R(0.  9)  =  1  + 

12.4  9  for  0  <  9  <  1°, 

(I-26b) 

and 

R( 0.  1)  -  R( 0. 9 )  = 

13.4  for  9  >  1°. 

(I-26c) 

The  distribution  functions  V  (p,  9)  and  V1  (p,  9)  are  related  by 

F  F 

VF(p,  9)  =  V'F(p,  9)  -  V'F(50,  9)  (1-27) 

Only  V  (p,  9),  which  has  the  median  of  V*  (p,  9)  as  a  reference 
r  F 

value,  was  used  in  determining  short-term  variability. 

It  was  unnecessary  to  consider  a  different  aircraft  antenna  gain 
for  the  direct  and  reflected  rays  because  their  arrival  angles  are 
nearly  equal  as  would  be  more  apparent  if  Fig.  1-1  were  drawn  to 
scale.  This  made  calculations  considerably  simpler  since  the  air¬ 
craft  antenna  gain  could  be  described  statistically  by  a  cumulative 

distribution  V  (p,  y  ),  which  is  a  function  of  y  only. 

.A.  OZ  OZ 

V  (p,  y  distributions  were  developed  using  data  from  an  air¬ 
craft  antenna  modeling  study,  employing  a  passenger  type  aircraft 
with  an  E-cavity  VOR  antenna  in  its  vertical  stabilizer  [Convair,  1959]  . 
The  data  used  were  gain  measurements  taken  around  the  model  with 

y  _  constant.  For  each  y  „  from  0°  through  -60°  the  measured 
o2  o2 

data  were  expressed  as  a  cumulative  distribution  and  this  distribution 
approximated  by  using  segments  of  normal  distributions  drawn  through 


I- 1 1 


the  1,  10,  50,  90,  and  99%  points  as  shown  in  Fig.  2  of  the  main  body 
of  the  report.  Linear  interpolation  was  used  to  obtain  distributions 
for  angles  at  which  data  were  not  available. 

Distributions  of  transmission  loss  values  were  calculated  as 
follows  using  the  terms  discussed: 

L(p,  d)  =  L  -  [V(p,  0  or  d)  *  V  (p,  0)  *  V  (p,  y  )]  (1-28) 

m  r  o  l. 

In  this  equation  all  terms  except  L  are  cumulative  distributions  for 

m 

the  reliability  p,  so  that  the  calculation  of  L(p,  d)  required  several 
successive  convolution  processes  which  are  indicated  by  the  operational 
symbol  and  will  be  further  discusses  later  on. 

2.  TACAN  PROPAGATION  MODEL 

The  TACAN  propagation  model  differs  from  the  VOR  model  dis¬ 
cussed  above  principally  by  the  use  of  vertical  instead  of  horizontal 
polarization  and  by  the  assumption  of  a  uniformly  rough  earth  model 
for  transmission  loss  calculations  at  the  1150  Mc/s  center  frequency 
of  the  TACAN  band.  The  vertical  radiation  pattern  of  the  ground 
station  was  shown  on  Fig.  3  of  the  main  body  of  this  report.  It  is 
based  on  a  typical  center  array  [Casabona,  19  56]  ,  with  the  antenna 
located  18  ft  above  the  center  of  the  VOR  counterpoise  and  30  ft  above 
ground. 

Median  transmission  loss  L  ,  again  including  the  gain  of  the 

transmitting  antenna,  was  calculated  using  geometric  optics  methods 

for  distances  within  the  radio  horizon,  but  the  reflected  energy  was 

treated  statistically.  Thus  L  becomes  mainly  a  function  of  the 

m 

direct  ray  path  distance  and  the  gain  of  the  transmitting  antenna  in 
the  appropriate  directions. 


1-12 


The  statistical  treatment  of  reflections  from  the  counterpoise 
was  based  on  the  assumption  that  the  complex  lobe  structure  due  to  the 
high  frequency  and  the  high  antenna  can  best  be  described  by  a  cumula¬ 
tive  distribution  function  obtained  by  adding  vectorially  two  rays,  one 
at  random  phase  relatively  to  the  other,  where  the  "constant"  vector 
represents  the  direct  ray,  and  the  vector  at  random  phase  has  a 
magnitude  r  times  the  magnitude  of  the  "constant"  vector  which  is 
given  by  the  assumed  reflection  coefficient  and  the  voltage  gains  of 
the  transmitting  and  the  receiving  antenna  as  a  function  of  the  grazing 
angle  and  the  other  auxiliary  angles  described  by  the  geometry  of  Fig. 
1-1.  This  distribution  function  can  be  considered  another  factor  enter¬ 
ing  the  short-term  variability,  as  the  movement  of  the  aircraft  through 
space  can  be  assumed  to  be  random  with  respect  to  the  complex  lobe 
pattern  described  by  the  distribution.  The  median  value  of  this  distri¬ 
bution  represents  a  term  added  to  the  free -space  median  L  calculated 

m 

as  a  function  of  ray  path  length,  and  the  distribution  itself  with  respect 
to  the  median  will  be  used  later  as  a  component  of  the  total  variability 
of  the  signal. 

This  concept  is  illustrated  on  Fig.  1-2.  The  cumulative  distribu¬ 
tion,  designated  V*  (p,  r),  due  to  the  vector  addition  of  a  unity  com- 
ponent  and  a  component  r  at  random  phase,  in  decibels,  is  given  by 
Norton,  Vogler,  Mansfield,  and  Short  [1955]: 

V'  (p,  r)  =  -10  log  [1  +  r2  +  2r  cos  t t  ( p/ 100)]  (1-29) 

where  p  is  the  percentage  of  time,  and  the  argument  tt  (p/ 100)  is  in 
radians.  The  median  of  this  distribution  is  obtained  from  setting  p 
equal  to  50: 
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V  q(  50,  r )  =  - 10  log  ( 1  +  r2)  (1-30) 

This  median  value  is  added  to  the  median  transmission  loss  L  , 

m 

and  takes  into  account  the  median  power  contained  in  the  ray  reflected 
from  the  counterpoise. 

The  short-term  variability  V  (p,  r)  due  to  this  distribution,  to 
be  used  in  the  same  manner  as  previously  shown  for  VOR,  is  given  by: 

VC(P»  r)  =  V'c(p,  r)  -  V'c(50,  r)  (1-31) 

The  relative  magnitude  r  of  the  random  component,  assumed 
to  be  0.9  for  the  illustrative  example  of  Fig.  I-Z,  is  given  by: 

r  =  (RCgrl)/gol  (I‘32) 

where  R  is  the  reflection  coefficient  of  the  counterpoise,  and  g 

O  1 

and  g  i  are  the  voltage  gain  factors  of  the  transmitting  antenna  rela¬ 
tive  to  an  isotropic  radiator  in  the  direction  of  the  direct,  and  of  the 
ground- reflected  ray,  respectively;  similar  to  the  VOR  case.  The 
dependence  of  the  voltage  gain  on  the  elevation  angles  y  ^  and  y  ^ 
(defined  previously  for  the  VOR  case  and  on  Fig.  1-1)  is  shown  for  the 
TACAN  ground  antenna  on  Fig.  3  of  the  main  body  of  the  report. 

Values  for  the  reflection  coefficient  R  were  assumed  as 
follows  as  functions  of  the  grazing  angle  ip : 

For  ip  <30°,  R  =  0  (no  reflection  from  the  counterpoise) 

Li 

For  30°  <  ip  <  43°,  R  =4.19  4'“  2.  19,  where  p  i-s  i n  radians 

L/ 

For  i|i  >43°,  R  =  0.  95 

V-> 
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Reflection  from  the  uniformly  rough  ground  was  assumed  to 
occur  for  grazing  angles  less  than  30  degrees.  In  this  case  a  cumula¬ 
tive  distribution  function  V'  (p,  K)  represents  the  energy  scattered 

F 

from  the  ground,  similar  to  the  function  V'  (p,  r)  used  to  represent 

reflections  from  the  counterpoise.  Here,  the  parameter  K  =  ZOlog^k 

where  k  is  the  ratio  of  the  root-mean- square  value  of  the  field  due  to 

scattering  of  the  reflected  ray  to  that  due  to  the  direct  ray.  With  an 

assumed  effective  ground  reflection  coefficient  R  ,  and  antenna 

G 

voltage  gain  values  g  and  g  ^  as  before,  k  is  given  by: 


k  =  (RGgrl,/gol  (I~ 

The  distribution  function  V'  (p,  K)  was  discussed  and  tabulated  by 

F 

Norton,  Vogler,  Mansfield,  and  Short  [1955]  for  various  values  of  K. 
The  effective  ground  reflection  coefficient  R^.  was  assumed  as 
follows  as  a  function  of  the  grazing  angle  ijj : 


For 
F  or 
For 


ij;  <30°,  R  = 

u 

30°  <  lJj  <  43°, 
4;  >  43°,  Rg  = 


0.  3 

R  =  -0.  13Z  lJj  +  0.  99Z,  where  i|i  is  3n  radians 
G 

0  (no  reflection  from  the  ground  due  to  the 
presence  of  the  VOR  counterpoise) 


From  the  distribution  of  V'  (p,  K),  the  median  value  for  p  =  50  was 

r 

determined  and  subtracted  from  the  median  transmission  loss  value 

L  in  order  to  take  into  account  the  median  power  reflected  by  the 
m 

rough  earth.  The  short-term  variability  V  (p,  K)  due  to  this  distri¬ 
bution,  and  to  be  used  later  on,  is  given  by: 


VF(p,  K)  =  V'F(P,  K)  -  V'f(50,K) 


(1-34) 
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The  geometric  parameters  according  to  Fig.  I- 1  are  the  same 
for  TACA.N  as  for  the  VOR  case  discussed  above.  The  treatment  of 
reflections  from  the  ground  eliminates  the  need  for  calculating  a  diver¬ 
gence  factor.  For  reflection  from  the  counterpoise  the  divergence 
factor  is  assumed  to  be  unity,  as  before.  The  carrier  frequency  used 
in  the  calculations  was  1150  Me/ s,  and  the  height  of  the  ground  antenna 
is  30  ft  above  ground  and  18  ft  above  the  counterpoise.  The  maximum 
gain  of  the  ground  antenna  is  8.  15  db  relative  to  an  isotropic  antenna. 
This  value  occurs  at  an  elevation  angle  of  approximately  5  degrees 
(see  Fig.  3  of  the  main  body  of  the  report,  where  the  ordinate  scale 
designates  the  voltage  gain  g  ^  relative  to  this  maximum). 

The  median  transmission  loss,  L  ,  for  the  TACAN  case,  in 

m 

view  of  the  procedures  described  above,  is  now  given  by: 

Lm  =  36.  58  +  20  log10  fMc  +  20  log10  ^  -  20  log1()  gQl 

-8.15+  V'  (50,  r)  -  V'  (50,  K)  (1-35) 

C  F 

Corrections  for  excess  ray  bending  were  made  in  the  same  way 
as  for  the  VOR  model,  using  the  same  exponential  atmosphere. 

Calculations  of  median  transmission  loss  beyond  the  radio 
horizon  were  made  by  combining  smooth-earth  diffraction  with  tropo¬ 
spheric  scatter  in  the  same  way  as  was  done  for  VOR.  However,  the 
method  of  calculating  diffraction  according  to  Vogler  [1961]  could  not 
be  used,  because  it  does  not  apply  to  vertical  polarization.  It  was 
replaced  by  more  general  methods  described  by  Norton,  Rice,  and 
Vogler  [1955],  and  in  NBS  Report  6767  [1961]  .  Again,  for  beyond- 
the -horizon  calculations  the  grazing  angles  are  essentially  zero,  and 
the  effect  of  the  counterpoise  is  neglected.  The  construction  of  final 
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propagation  curves  (L  versus  d)  for  various  discrete  aircraft  alti¬ 
tudes  involved  again  a  certain  amount  of  blending  from  one  propagation 
mechanism  to  another. 

Long-term  variability  was  described  statistically  for  all  distances 
by  the  function  V(p,  d),  using  the  same  methods  as  in  the  VOR  propa¬ 
gation  model.  For  the  TACAN  model  the  function  V(p,  9)  was  not  used 
at  all 

Short-term  ( within-the  -hour )  fading  within  the  radio  horizon  was 
described  by  cumulative  distribution  functions  V  (p,  r)  and  V  (p,  K), 
which  result  from  the  random  addition  of  the  direct  and  the  counter¬ 
poise  or  ground-reflected  ray  as  described  above.  Equations  (1-31) 
and  (1-34)  defined  these  functions.  For  distances  beyond  the  radio 
horizon,  V  (p,  K)  was  used  with  K  determined  in  a  way  similar  to 
the  one  described  for  the  VOR  model,  but  based  on  1,  000  Me/ s  fading 
range  data  from  the  paper  by  Janes  [1955]  .  The  two  curves  for  K 
obtained  were  blended  in  the  vicinity  of  the  radio  horizon  in  order  to 
obtain  a  continuous  function  of  K  versus  distance. 

As  in  the  case  of  VOR,  variations  in  the  gain  of  the  aircraft 
antenna  also  enter  as  a  short-term  fading  term,  and  are  described  by 
a  cumulative  distribution  function  V  (p,  y  ),  where  y  is  the 

■/a.  O  i—i  O 

arrival  angle  of  the  direct  ray  at  the  aircraft  (see  Fig.  1-1).  For  the 
type  of  aircraft  antenna  used  in  these  calculations,  typical  distribu¬ 
tions  for  various  values  of  y  are  shown  on  Fig.  4  of  the  main  body 

o2 

of  the  report. 

All  time  variability  functions  were  combined  by  statistical  con¬ 
volution,  and  the  resulting  cumulative  distribution  of  transmission 
loss  as  a  function  of  the  time  variability  p  and  the  distance  d  is 
given  by: 
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L(p,ci)  =  Lm  -  {y(P’ d)  *  VA(p,  Vo2)  *  VF(P>  K)  *  [-Vc(p,  r )]  }  (1-36) 

All  terms  in  (1-36)  except  L  are  cumulative  distributions,  so  that 

m 

several  successive  convolutions  are  required,  which  are  denoted  by 
the  symbol  *.  The  choice  of  the  values  for  the  reflection  coefficient 
from  the  ground  and  the  counterpoise  assured  proper  use  of  the  com¬ 
ponents  V  (p,  K)  and  V  (p,  r):  for  grazing  angles  less  than  30  degrees 
the  cumulative  distribution  V  (p,  r)  degenerates  into  a  constant. 

3.  COMPUTATION  PROCEDURES 
a.  General  Methods  and  VOR  Calculations 

The  model  on  which  the  computation  methods  and  the  numerical 
results  in  this  report  are  based  considers  levels  of  signals  from  the 
desired  and  from  the  interfering  station,  with  the  receiving  aircraft 
located  directly  above  the  great  circle  path  between  the  two.  First, 
the  computer  program  applicable  to  the  VOR  system  will  be  discussed 
in  detail  as  representative  of  the  overall  approach.  Variations  intro¬ 
duced  by  consideration  of  the  TACAN  system  will  be  pointed  out  later. 
Although  the  program  has  been  written  in  a  way  so  that  it  is  applicable 
in  the  future  to  a  multi- station  interference  study,  its  organization  is 
better  understood  if  the  discussion  is  limited  first  to  the  two-station 
model  described  above.  Once  this  has  been  understood  possible  exten¬ 
sions  of  the  program  to  more  complicated  situations  are  easier  to 
follow. 

The  machine  program  consists  essentially  of  three  parts  which 
may  be  taken  as  three  different  "levels"  of  control.  Actual  computa¬ 
tions  are  performed  in  Levels  I  and  II,  with  Level  III  being  a  catalogue 
of  computer  routines  available  on  call.  Level  II  acts  principally  to 
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insure  proper  sequencing  of  the  routines  involved  and  keeps  track  of 
the  operations  performed.  Fig.  1-3  is  a  block  diagram  showing  the 
relations  between  these  levels  and  the  sequence  of  computations. 

Level  I  uses  as  its  input  a  set  of  basic  physical  data  in  the  form 
of  tables,  which  have  been  previously  assembled  by  separate  machine 
programs  and  other  computations.  These  basic  data  include  transmis¬ 
sion  loss  as  a  function  of  great-circle  distance  and  aircraft  height  in 
the  way  defined  in  the  previous  sections,  together  with  other  geometric 
and  equipment  parameters  necessary  for  the  problem.  The  spatial 
relationship  of  the  ground  stations  involved  is  also  part  of  the  basic 
data,  but  the  aircraft  antenna  pattern  (given  by  its  cumulative  gain 
distribution  V(p,  y  ,)  is  brought  in  only  at  Level  II. 

It  has  been  noted  before  that  in  this  application  the  term  trans¬ 
mission  loss  includes  the  gain  of  the  transmitting  ground  antenna  as 
well  as  the  effect  of  reflections  from  the  ground  or  the  counterpoise. 

The  construction  of  the  basic  data  tables  is  accomplished  independently 
of  the  main  computer  program  and  involves  a  combination  of  graphical 
and  computer  methods.  The  tables  contain  as  the  argument  the  distance 
d  along  the  great  circle  path  on  the  ground  from  the  transmitting 
station  considered  to  the  aircraft  location,  and  as  functional  values 
various  geometric  parameters  defined  in  Fig.  I—  1 ,  the  angular  dis¬ 
tance  9  and  other  parameters  used  in  the  calculation  of  the  median 

transmission  loss,  L  .  A  separate  table  is  used  for  each  aircraft 

m 

height  considered;  the  specific  heights  used  were  listed  in  Section  6 
of  the  main  body  of  this  report.  Although  the  formulation  of  these 
tables  is  an  involved  process,  they  are  useful  as  they  ultimately  result 
in  a  saving  of  machine  time  in  the  main  program.  Once  a  table  of 
transmission  loss  values  is  available,  the  machine  calculation  of  a 


1-19 


value  for  a  specific  distance  is  accomplished  by  a  table  look-up  and  an 
interpolation. 

The  characteristics  of  the  tables  dictate  the  following  order  of 

calculations  in  Level  I,  as  illustrated  on  Fig.  1-4.  Tables  of  all  basic 

data  for  a  specific  value  of  aircraft  altitude  are  read  into  storage.  Then 

input  information  is  read  which  specifies  the  separation  S  =  (d^  +  d^), 

where  d^  and  d  ,  respectively,  are  the  distances  along  the  great 

circle  path  on  the  ground  from  the  desired  and  the  undesired  station  to 

the  aircraft.  An  initial  value  for  d  is  contained  in  the  input  and 

appropriate  information  regarding  increments  is  part  of  the  program. 

In  Level  I,  the  median  transmission  loss  value  L  and  other  para- 

m 

meters  are  calculated  for  d  and  d  ,  and  presented  to  Level  II  for 
further  calculations.  Thus  control  alternates  between  Level  I  and 
Level  II  for  each  aircraft  position  with  respect  to  the  two  ground 
stations  at  the  given  specific  altitude.  The  schedule  is  repeated  for 
all  desired  values  of  transmitter  separation  at  this  altitude;  subse¬ 
quently  new  tables  for  a  different  altitude  are  read  into  the  machine, 
and  the  entire  procedure  repeated  until  all  available  input  values  are 
exhausted. 

The  parameters  calculated  by  Level  I  of  the  program  are  used 
in  Level  II  for  the  calculation  of  the  various  cumulative  distributions 
which  make  up  the  expression  for  the  time  distribution  of  transmission 
loss  previously  shown  as  (I-Z8): 

L(p,  d)  =  L  -  [V(p,  0  or  d)  *  V  (p,  0)  *  V  (p,  y  )]  (1-28) 

m  r  o  c 

Here  the  term  V(p,  0  or  d)  is  the  cumulative  distribution  repre¬ 
senting  long-term  variability  of  hourly  medians  (in  the  VOR  case, 
either  as  a  function  of  the  angular  distance,  0,  for  beyond-the-horizon 
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paths,  or  as  a  function  of  the  path  distance,  d,  for  within-the -horizon 
paths).  V  (p,  9)  represents  the  short-term  (within-the -hour )  variability 
and  V  (p,  y  7)  represents  the  variability  introduced  by  the  aircraft 
antenna  as  a  function  of  the  elevation  angle  y  defined  on  Fig.  1-1. 

The  percentage  of  time  is  represented  by  the  symbol  p,  and  L(p,  d) 
denotes  the  resulting  cumulative  distribution  of  transmission  loss.  In 
this  program,  L(p,  d)  is  calculated  for  both  distances  d  and  d^ 
defined  above . 

The  symbol  *  in  (1-28)  denotes  the  statistical  convolution  of 

cumulative  distributions.  Basically,  if  two  variables  x  and  y  are 

statistically  independent  and  given  by  their  cumulative  distributions, 

their  statistical  convolution  is  the  cumulative  distribution  of  either  the 

variable  z  =  x  +  y,  or  the  variable  z1  =  x  -  y,  depending  on  the 

application  of  the  result,  and  determined  in  the  machine  program  by 

an  appropriate  input  code.  The  cumulative  distribution  of  z  is 

obtained  most  simply  by  selecting  n  equally  spaced  percentage  values 

from  the  individual  distributions  of  both  x  and  y  calculating  all 

possible  sums  z,  =  x.  +  y.,  and  forming  the  cumulative  distribution 
k  i  j 

of  all  values  of  z  obtained  in  this  manner.  Similarly,  the  cumula- 

K 

tive  distribution  of  z'  is  obtained  by  forming  all  possible  differences 

z.  '  =  x.  -  y.. 
k  1  J 

For  the  convolution  of  the  individual  distributions  used  in  this 
study,  a  somewhat  more  sophisticated  method  was  used  in  selecting 
percentage  intervals,  which  provides  greater  accuracy  at  the  "tails" 
of  the  distributions.  All  individual  distributions  are  contained  in 
storage  in  tabular  form,  and  the  convolution  process  itself  takes  the 
form  of  a  program  routine.  Successive  convolutions  (as  in  1-28)  are 
handled  by  first  convoluting  two  distributions  and  then  using  this  inter¬ 
mediate  result  in  a  convolution  with  the  third  individual  distribution. 
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The  convolution  process  is  commutative;  i.e.,  the  order  in  which  the 
individual  distributions  are  used  does  not  matter. 

For  each  given  aircraft  elevation  and  spacing  of  ground  station, 
the  distribution  of  transmission  loss  L(p,  d)  is  determined  for  d  =  d 
and  d  =  d^,  the  distances  of  the  aircraft  from  the  two  ground  stations 
considered.  From  the  distribution  L(p,  d^),  which  is  determined  first, 
the  values  L(0.05,  d^)  and  L(0.95,  d  )  are  extracted,  converted  to 
relative  voltages  E(0.05,  d^)  and  E(0.95,d^)  by  a  special  program 
routine,  and  stored  for  later  use.  Next,  the  transmission  loss  distri¬ 
bution  L(p,  d  )  corresponding  to  the  path  from  the  desired  station  to 
the  aircraft  is  determined.  Finally,  the  cumulative  distribution  of  the 
difference  in  transmission  loss  values,  L(p,  d^.)  *  [-L(p,  d  )]  is 
obtained  by  a  convolution.  As  discussed  above,  this  is  equivalent  to 
the  distribution  of  the  desired  to  the  undesired  signal  strength  ratios 
D/U  at  the  receiver  input  in  the  aircraft.  From  this  last  distribution, 
a  number  of  percentage  values  between  D/U(0.01)  and  D/U(0.99) 
are  selected  and  stored.  They  are  also  punched  on  cards,  together 
with  the  values  d  and  d^  as  well  as  the  voltages  E(Q.  05,  d  )  and 
E(0.95,  dy),  previously  obtained,  which  designate  the  extent  of  the 
noise -limited  service  as  a  function  of  distance  and  aircraft  height. 

The  subscript  U  in  this  case  has  no  further  significance;  it  only 
denotes  that  the  data  for  the  path  identified  by  U  were  used  for  this 
purpose.  The  output  in  the  form  of  punched  cards  facilitates  graphical 
presentation  of  the  results  using  an  automatic  plotter  with  card  input. 

Many  of  the  distribution  functions  which  enter  into  the  machine 
calculation  of  L(p,  d)  are  empirical  and  can  be  stored  in  the  machine 
in  tabular  form  only.  For  the  purpose  of  consistency  and  simplicity, 
all  distribution  functions  entering  the  program  routines  are  handled  in 
tabular  form,  even  if  they  can  be  expressed  analytically.  Consequently 
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a  large  number  of  tables  is  involved  and  a  considerable  amount  of 
"bookkeeping"  is  necessary,  which,  together  with  the  sequencing  of 
computation  and  output  routines,  is  the  principal  function  of  the  Level 
II  control. 

Fig.  1-5  is  a  flow  diagram  of  Level  II  computations  and  illustrates 
the  interplay  between  Level  II  and  the  program  routines  contained  in 
Level  III.  These  routines  are  designated  by  code  names  on  Fig.  1-5, 
which  are  tabulated  and  further  explained  below. 


Table  of  Program  Routines  Contained  in  Level  III 


Code  Designation 


Function 


(a)  Major  Routines 


ZEBRA 


Computes  a  table  describing  the  cumulative 
distribution  function  of  the  receiving  antenna 


pattern,  VA(p,  yo2). 


FALTG 


Calculates  the  convolution  of  cumulative 
distribution  functions;  results  are  available 
to  the  program  in  form  of  a  table. 


VPT 


Calculates  the  cumulative  distribution 
function  V(p,  9);  results  are  in  tabular 
form. 


VPD 


Calculates  the  average  cumulative  distri¬ 
bution  function  V(p,  d)  from  summer  and 
winter  values  available  from  BRAP;  results 
are  in  tabular  form. 
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Table  of  Program  Routines  Contained  in  Level  III  (Continued) 


Code  Designation 


(b)  Subroutines 
INTERP 


ENDPNT 


Function 


This  is  an  interpolation  program,  using  the 
n^1  Lagrangian  approximation,  with  n 
specified  in  any  particular  case.  It  is  used 
as  a  part  of  ZEBRA  and  FALTG,  where 
interpolation  between  tabulated  values  is 
necessary. 

This  is  a  routine  similar  to  an  interpolation, 
but  it  fits  Gaussian  distribution  functions 
between  given  points  of  cumulative  distri¬ 
butions  like  V(p,  d),  etc.  It  is  used  in 
ZEBRA  and  VPD. 


BRAP  This  routine  handles  the  actual  calculation 

of  summer  and  winter  values  for  V(p,  d) 
at  selected  values  of  p.  It  is  used  for  VPD, 
but  ENDPNT  is  needed  to  find  intermediate 
value  s. 


Summarizing  the  program  steps,  it  is  seen  that  the  group  of  pro¬ 
grams  contained  in  Level  I  are  called  by  Level  II  to  carry  out  the  actual 
computations  of  transmission  loss  values  and  necessary  geometric 
parameters.  Level  II  control  makes  decisions  based  on  information 
furnished  by  Level  I  as  to  which  programs  of  that  level  are  pertinent; 
subsequently  Level  II  control  calls  these  programs  in  their  proper 
sequence,  provides  them  with  proper  input,  and  keeps  track  of  their 
output.  This  output  eventually  becomes  input  for  additional  Level  III 
program  steps  also  controlled  by  Level  II,  unless  it  constitutes  and 
end  result  to  be  read  out  of  the  machine. 
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There  is  a  specific  routine  for  each  component  distribution  which 
can  enter  into  the  transmission  loss  distribution  given  by  (1-28).  Upon 
being  given  a  particular  value  of  the  elevation  angle  y  the  program 
for  the  antenna  gain  distribution  function  VA(p,  v  )  extracts  values  of 
this  function  for  p  =  1 ,  10,  50,  90,  and  99%  from  the  tables  contained 
in  storage.  These  and  the  corresponding  percentage  values  used  are 
presented  to  the  ENDPNT  routine  which  constructs  a  new  table  by  a 
process  equivalent  to  plotting  the  available  points  on  normal  probability 
paper  and  drawing  straight  lines  between  them.  Thus  a  more  complete 
table  of  the  distribution  function  is  obtained  for  use  in  the  calculations. 

The  calculation  of  the  long-term  variability  function  V(p,  d)  uses 
a  similar  procedure,  but  in  this  case  two  distributions  V^(p,  d)  and 
V^p,  d)  must  first  be  calculated,  corresponding  to  summer  and  winter 
data  [NBS  Report  6767]  .  These  two  distributions  are  averaged  to  give 
the  desired  V(p,  d)  by  means  of  the  BRAP  routine. 

Tables  are  used  exclusively  for  the  long-  and  short-term  variabi¬ 
lity  functions  V(p,  9),  and  V  (p,  0)  used  in  the  VOR  study. 

With  the  extensive  use  of  tables,  the  need  to  interpolate  arises 
often,  and  linear  interpolation  will  not  always  suffice.  Therefore,  a 
subroutine  was  included  in  the  program  which  calculates  the  n^1 
Lagrangian  approximation  [Lance,  I960],  where  n  :5~  m  -  1  with  m 
being  the  number  of  entries  in  the  table.  The  special  case  n  =  1  results 
in  linear  interpolation.  As  a  practical  matter,  n  is  usually  required 
to  be  less  than  five.  This  is  necessitated  by  the  fact  that  Lagrangian 
interpolation  may  exhibit  large  oscillations  between  tabular  entries 
whenever  they  do  not  lie  on  rather  smooth  analytic  functions. 
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b.  TACAN  Calculations 

The  program  for  the  TACAN  system  is  only  slightly  different 
from  the  one  discussed  above  for  VOR.  In  Level  I,  the  tables  contain¬ 
ing  data  for  use  in  transmission  loss  calculations  are  different  due  to 
the  changes  in  the  propagation  model  discussed  in  Section  2  of  this 
Appendix.  These  tables  contain  the  great-circle  distance  d  as  the 
argument  as  before,  but  in  addition  to  the  TACAN  median  transmission 
loss  L  the  parameters  r  and  K  defined  by  (1-32)  and  (1-33), 
respectively,  are  also  listed  as  functions  of  d.  These  parameters 
enter  two  component  distributions  V  (p,  K)  and  V  (p,  r)  used  for 
TACAN  which  represent  the  variability  introduced  by  the  scattered 
ground-  (or  counterpoise-)  reflected  power  and  by  the  multiple  lobes 

of  the  ground  antenna.  The  distribution  function  V  (p,  K)  replaces 

F 

V  (p,  9)  used  for  VOR.  The  long-term  variability  in  the  TACAN  case 
F 

is  always  represented  by  the  V(p,  d)  functions.  The  resulting  trans¬ 
mission  loss  distribution  for  the  TACAN  case  is  now  given  by  (1-36), 
which  contains  the  additional  term  -V  (p,  r): 

L(p,  d)  =  L(m)  -  {V(p,  d)  *  VA(p,  yo2)  *  VF(p,  K)  *  [-^(p,  r )]  }  (1-36) 

Level  II  differs  from  VOR  only  by  the  substitution  of  different  tables 
corresponding  to  the  TACAN  distribution  functions  mentioned,  and 
Level  III  is  changed  accordingly  by  deleting  routines  not  needed  and 
including  new  routines  for  the  calculation  of  the  distribution  functions 

V  (p,  r)  and  V  (p,  K).  Also,  the  tables  used  to  calculate  the  antenna 
gain  distribution  function  VA(p,  y  )  are  changed  in  accordance  with 
the  TACAN  aircraft  antenna  model  used.  For  TACAN,  the  noise- 
limited  service  (without  regard  to  interference  from  another  station) 
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is  denoted  by  receiver  input  power  levels  rather  than  by  voltage  levels 
as  used  for  VOR.  Consequently,  power  levels  P  (0.95,  d  )  and 

3.  U 

P(0.  05,  d  )  are  calculated  from  the  corresponding  transmission  loss 
values  and  stored  for  output. 

c.  Extension  of  the  Program  to  the  Case  of 
Multiple -Station  Interference 

The  present  study  is  limited  to  a  single  interfering  station  acting 
on  the  reception  of  the  desired  signal  at  the  aircraft  receiver  terminals. 
The  computer  program,  however,  can  be  expanded  to  permit  a  gener¬ 
alization  of  the  problem;  i.e.,  the  consideration  of  interference  from 
more  than  one  station.  In  this  case  it  would  be  necessary  to  provide 
more  input  information  than  given  by  the  station  separations  in  order 
to  specify  completely  their  positions  relative  to  each  other  and  relative 
to  the  aircraft.  This  requirement  affects  Level  I  routines  by  demand¬ 
ing  more  input  parameters  in  several  sets.  Furthermore,  an  additional 
input  to  Level  II  is  needed  which  indicates  the  number  of  interfering 
stations  involved. 

When  considering  n  stations,  Level  II  of  the  program  would 

first  calculate  L(p,  d^^),  L(P’du(Z))’ . Mp»  dU(n  -  1  )^s  which 

are  the  transmission  loss  distributions  corresponding  to  each  inter¬ 
fering  station.  These  must  then  be  converted  to  relative  power  units, 
and  the  distribution  of  their  sum  is  determined  by  successive  convolu¬ 
tions  under  the  simplified  assumption  that  no  correlation  exists  between 
the  various  distributions.  The  result  is  converted  back  to  a  transmis¬ 
sion  loss  distribution  and  substituted  for  the  L(p,  d^)  distribution  used 
in  the  single  -  inte  rfe  rence  case.  No  other  changes  are  necessary  in 
Level  II,  neither  is  Level  III  affected  by  the  number  of  stations  considered. 


GEOMETRY  FOR  RAY  INTERFERENCE  WITHIN  THE  RADIO  HORIZON 
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THE  CUMULATIVE  DISTRIBUTION  FUNCTION  v'c(p,  r)  FOR  r  =  0.9 
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BLOCK  DIAGRAM  OF  COMPUTATION  PROCEDURE 
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Figure  1-3 
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FLOW  DIAGRAM  FOR  LEVEL  I  OF  MAIN  COMPUTER  PROGRAM 


rom  Level  II  To  Level  II 
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Figure  1-5 
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